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1.  PREFACE  AND  TIMETABLE 


This  report  presents  progress  during  the  first  year,  and 
outlines  modifications  in  the  research  plan  of  the  second  year  of  a 
three-year  continuation  of  our  ongoing  research  effort.  This  research 
seeks  to  predict  shelf  and  beach  forms  and  their  changes  from  knowledge 
of  local  bathymetry  and  the  driving  forces  due  to  winds,  waves  and 
currents  and  their  complicated  interactions  with  nearshore  sediments. 

The  full  three-year  effort  was  detailed  in  the  proposal  submitted  last 
year. 

Table  1  shows  the  research  timetable  proposed.  There  are  no 
major  changes  in  the  planning.  The  following  progress  reports  cover 
the  individual  tasks  addressed  in  the  1980-81  contract  year,  as  indicated 
by  arrows. 


2.  GENERAL  RESEARCH  PLAN  1980  -  1983 

The  proposed  work  falls  into  three  distinct  but  interdependent 
areas  of  research.  They  are: 

A.  Wave  and  current  dynamics, 

B.  Fluid-sediment  interactions, 

C.  Coastal  zone  remote  sensing. 

Each  area  logically  divides  into  several  "tasks"  addressing  specific 
scientific  or  instrumentation  problems.  These  tasks  are  also  interrelated, 
but  vary  in  complexity  and  duration  depending  on  the  objectives.  They  will 
be  phased  in  and  out  of  the  research  program  depending  on  the  scientific 
merit  and  results,  the  interest  of  ONR  and  the  interests  and  goals  of 
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the  principal  investigators.  The  proposed  tasks  and  their  tentative 
timetables  for  the  three  year  period  1980-1983  are  shown  in  Table  1. 

It  is  clear  that  the  overall  objective  of  the  research  program 
dictates  the  specific  goals  of  the  individual  tasks.  For  example,  waves 
and  currents  in  nearshore  waters  are  the  driving  forces  for  sedi¬ 
ment  motion,  as  well  as  being  interesting  and  important  in  their 
own  right.  Progress  on  the  prediction  of  shelf  and  beach  topogra¬ 
phic  forms  clearly  requires  improvements  in  the  specification  of 
fluid-driving  forces.  Fluid-sediment  interactions  are  the  corner¬ 
stone  of  this  research  project,  and  also  pose  the  most  difficult 
and  complex  set  of  problems.  Even  the  most  fundamental  laws,  the 
equations  of  motion,  for  granular  fluids  are  not  completely  known, 
much  less  solved.  The  analogous  laws  for  fluids  have  been  known 
for  hundreds  of  years.  Approximations  and  empirical  deductions, 
however,  still  allow  progress  to  be  made  in  many  areas  of  fluid- 
sediment  interaction.  Frequently,  fluid  and  fluid-sediment  inter¬ 
action  experiments  are  run  simultaneously,  emphasizing  the  synergy 
between  these  objectives. 

The  objective  of  the  remote  sensing  task  is  to  use  and  refine 
existing  techniques  to  study  beaches  and  shelves,  rather  than  to 
develop  new  remote  sensing  instrumentation.  *The  opportunities  for 
wide  spatial  coverage,  without  great  time  lags,  of  course  makes  re¬ 
mote  sensing  an  attractive  way  of  gathering  data  which  would  be  pro¬ 
hibitively  expensive  to  obtain  using  traditional  "ground  truth"  sen¬ 


sors. 


Table  1 

FLUID-SEDIMENT  INTERACTION  ON  BEACHES  AND  SHELVES 
Research  Timetable 
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980-81  1981-82  1982-83 


A.  WAVE  AND  CURRENT  DYNAMICS 
"AIR-SEA  DRIVING  FORCES" 

1.  Incident  Wave  Fields 

1.1  Spectral  Moments 

1.2  Island  Shadowing 

1.3  Shoaling  Waves 

2.  Shelf  Waves 

3.  Surf  Zone  Dynamics 

3.1  Turbulence 

4.  Current  Measurements  by  Geomagnetic  Induction 

5.  Other  Driving  Forces 

B.  FLUID-SEDIMENT  INTERACTIONS 
"BED  RESPONSE" 

1.  Sand  Bed  Response  to  Fluid  Stress 

1.1  "Oscillatory"Bursting 

1.2  Rheology  of  Granular  Fluid  Materials 

1.3  Bedload  Meter  Deployment 

2.  On-Offshore  Transport  of  Sand 

3.  Coastal  Forms 

3.1  Bedforms  Produced  by  Fluid-Granular  Interaction 

3.2  Cusps 

4.  Shore  Processes  of  Northern  Adriatic  Beaches 

C.  COASTAL  ZONE  REMOTE  SENSING 


1.  Remote  Sensing  of  Waves  (SAR) 

2.  Shelf  and  Nearshore  Dynamics 


The  underlying  assumption  behind  this  research  project  is  that 
the  dynamics  of  the  shore  and  shelf  areas  of  the  world  are  governed 
by  the  same  physical  processes,  although  individual  processes  may 
differ  in  importance  in  specific  areas  thus  producing  the  world's 

varied  coastlines.  Much  of  coastal  research  is  directed  toward  the 
detailed  description  of  specific  areas  and  the  treatment  of  these  areas 
as  unique  entities.  In  contrast,  this  work  seeks  to  understand  these 
basic  underlying  physical  processes  so  that  the  results  may  be  applied 
to  all  coastal  areas. 


A.  WAVE  AND  CURRENT  DYNAMICS 
"AIR-SEA  DRIVING  FORCES" 


The  principal  area  of  interest  in  the  broad  sense  continues  to 
be  wind  generated,  surface  gravity  waves  and  currents  and  other  fluid 
modes  which  derive  their  energy  from  these.  Emphasis  over  the  past  year 
has  been  on:  developing  estimation  techniques  for  momentum  flux 
(Task  Al.l);  completion  of  the  island  shadowing  project  (A1.2);  develop¬ 
ment  and  testing  of  a  broadly  applicable  wave  shoaling  model  (A1.3);  an 
analysis  of  run-up  and  swash  data  at  surf-beat  frequencies  (A2);  prelim¬ 
inary  work  to  measure  the  gross  properties  of  turbulence  in  and  near  the 
surf  zone  (A3.1);  and,  in  addition,  instrument  development  has  proceeded 
in  preparation  for  measurements  of  the  water  flow-induced  electric 
potential  across  the  entrance  of  a  tidal  basin  (A4). 

The  ultimate  interest  in  studying  the  forcing  functions  is  to 
understand  their  mechanics  so  that  their  full  potential  as  driving  forces 
may  be  applied  to  the  study  of  fluid-sediment  interactions  and  to  sediment 
response.  This  is  the  guiding  principle  in  selecting  areas  of  study. 


Al.  Incident  Wave  Fields 

Al.l  Spectral  Moments  (S.  S'.  Pawka) 

The  onshore  flux  of  longshore  directed  wave  momentum,  S  ,  is 

yx 

an  important  forcing  function  for  nearshore  processes.  S  has  been 

yx 

shown  to  drive  the  longshore  current  inside  the  surf  zone  (Bowen,  1969) 
and  to  be  an  Important  factor  in  the  longshore  transport  of  sand  (e.g. 
Inman  and  Bagnold,  1963;  Ballard  and  Inman,  1981).  Coastal  researchers 
have  expended  considerable  effort  in  attempts  to  make  accurate 


measurements  of  this  wave-field  parameter.  However,  the  statistical 
nature  of  these  measurements  has  not  been  previously  addressed. 


S  ,  expressed  as  a  function  of  frequency,  is  related  to  the 
yx 

following  moment  of  the  frequency-directional  spectrum,  E(f,a), 

2  * 

S  (f)  =  n(f) f  cosa  sina  E(f,a)  da 

y* 


(Al.l) 


where  a  is  propagation  direction,  f  is  the  wave  frequency,  and  n(f) 

is  the  ratio  of  group  to  phase  speed.  Accurate  estimates  of  S  (f)  can 

be  computed  from  accurate  measurements  of  E(f,a),  say  from  a  pressure 

sensor  array,  and  integrating  as  in  Al.l.  In  principle,  direct  estimates 

of  the  SyX  moment  can  be  made  with  a  "slope  array"  (Seymour  and 

Higgins,  1977)  or  a  biaxial  current  meter  (Pawka,  in  prep.).  For  example, 

the  co-spectrum  of  the  slopes,  n  =  I7  »  n ,  *  obtained  from  the  slope 

x  oX  y  ay 

array  has  the  form 


C_  n  (f)  *  k^  f  cosa  sina  E(f,a)  da  (A1.2) 

\ny  0 

These  direct  measurement  techniques  have  two  limitations.  First, 
these  "orthogonal  component"  systems  have  relatively  poor  resolution  power 
for  the  definition  of  the  directional  spectrum  E(f,a).  The  poor  reso¬ 
lution  would  render  these  systems  inadequate  for  the  detailed  study  of 
the  directional  wave  properties  discussed  in  Section  A1.2.  Second,  the 
(Gaussian)  statistical  nature  of  waves  leads  *to  a  variance  of  an  as  yet 

undetermined  form  in  the  estimates  of  S  .  These  statistical  variations 

yx 

and  their  effect  on  "S  array"  data  applications  are  addressed  in  this 

yx 

task. 


The  probability  density  function  for  the  co-spectrum  in  a 
2-dimensional  Gaussian  random  field  was  developed  by  Goodman  (1957).  Three 
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example  density  functions  of  this  distribution  family  are  shown  in 
Figure  Al-1.  The  effect  of  varying  directional  spectrum  characteris¬ 
tics  and  sampling  specifications  on  the  shape  of  these  distributions 
was  studied.  The  dependence  of  distribution  variance  and  skewness 
on  the  width  of  a  unimodal  directional  peak  and  the  degree  of  spectral 
averaging  is  shown  in  Figure  Al-2.  Both  parameters  strongly  influence 
these  low  order  distribution  moments.  Broad  directional  spectra  increase 

the  expected  variability  of  S  (f)  estimates  requiring  longer  sampling 

yx 

periods  for  stable  measurements.  Using  the  analytic  density  functions, 

appriori  knowledge  of  the  wave  directional  characteristics  can  be  used 

in  the  design  of  sampling  specifications  of  nearshore  experiments  or  wave 

climatological  measurement  systems. 

The  variance  of  S  (f)  estimates  is  lowest  for  unidirectional 
yx 

wave  conditions  (Figure  Al-3).  Estimates  of  zero  mean  S  (f)  conditions 

yx 

may  have  relatively  large  absolute  variances  while  this  is  not  true  for 

frequency  spectrum  estimates.  The  directional  moments  which  determine 

the  variance  of  S  (f)  estimates  can  all  be  directly  measured  with  a 
yx 

biaxial  current  meter  system.  Therefore,  this  limited  system  not  only 

yields  an  estimate  of  S  (f )  but  also  its  expected  variance. 

yx 

A  linear  array  of  sensors  can  be  used  to  estimate  the  wave 

directional  spectra.  The  integration  of  these  spectra  (Equation  Al.l) 

yields  an  estimate  of  S  (f).  The  probability  density  functions  discussed 

yx 

above  are  not  strictly  applicable  to  array  measurements,  which  have  some 
spatial  averaging  of  the  wave  field.  However,  Pawka  (in  prep.)  has  demon¬ 
strated  good  comparison  of  linear  array  estimates  of  S  (f)  with  those 

yx 

obtained  from  a  slope  array.  These  comparisons  included  close  agreement 
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Figure  Al-2.  Dependence  of  normalized  2nd  and  3rd  order  moments  of 

S  ( f )  distribution  on  the  width  of  a  unimodal 
yx 

directional  spectrum  and  the  number  of  ensembles 
averaged.  The  directional  distribution  employed 

was’  E(a)  =  cosm(a-10°)  where  o  is  the  propaga¬ 
tion  direction. 


VAR  (Syx(  f )) ;  DATA 
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VAR  (Syx  ( f ))  ;  THEORY 


Figure  Al-3.  Comparison  of  theoretical  values  of  the  variance 
of  SyX(f)  estimates  with  sample  values  obtained 

from  linear  array  data.  The  x's  represent  array 

data  which  is  normalized  to  the  energy  density  of 

a  single  sensor  while  the  dots  are  normalized  to 

an  average  of  3  sensors.  The  variances  have  the 

4  2 

units  cm  sec  . 


''A 


of  the  estimator  variance.  Therefore,  this  distribution  family  may  be 

useful  in  the  description  of  array  estimate  variability. 

A  400  m  length  linear  array,  described  by  Pawka  et  al  (1980), 

was  maintained  at  Torrey  Pines  Beach,  California.  Sample  distributions 

of  S  (f)  estimates  obtained  from  the  array  data  are  shown  in  Figure  Al-1. 
yx 

These  sample  distributions  are  compared  to  analytic  density  functions 
which  were  generated  using  the  measured  directional  spectra  from  the 
array.  The  comparisons  show  good  agreement  of  the  low  order  distribu¬ 
tion  moments.  The  analytic  and  sample  variances  (Figure  Al-3)  and  skew¬ 
ness  show  very  high  correlation.  The  comparison  of  higher  order  moments 
were  inconclusive  due  to  the  limited  nature  of  the  data  sets. 

Additional  field  data  for  comparisons  with  the  theoretical 
distributions  were  sampled  during  the  past  year.  Two  biaxial  electro- 
megnetic  current  meters,  separated  by  90  m,  were  placed  in  a  line  along 
the  5  m  depth  contour  at  the  Torrey  Pines  Beach  site.  The  400  m  array, 
located  in  a  mean  depth  of  9.6  m,  was  also  operated  for  background  wave 
information.  Roughly  5  days  of  long  data  runs  (duration  M5  hours)  were 

sampled.  This  data  will  be  used  to  further  test  the  S  (f)  distributions 

yx 

at  a  single  location  and  also  the  covariability  of  S  (f)  estimates  made 

yx 

at  a  given  spatial  separation. 
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A1.2  Island  Shadowing  (S.  S.  Pawka  and  D.  L.  Inman) 

Torrey  Pines  Beach,  Southern  California,  is  partially  sheltered 
from  the  deep  ocean  by  the  Channel  Islands.  Last  year's  progress  report 
described  field  data  comparisons  with  an  island  blocking  and  refraction 
model  of  the  wave  transformations  from  deep  ocean  to  the  coast.  The  islands 
were  shown  to  significantly  affect  the  coastal  wave  conditions.  However, 
only  limited  confirmation  of  the  blocking/refraction  model  in  the  mid  fre¬ 
quency  range  (0.082-0.114  Hz)  was  obtained  due  to  the  poor  quality  of  the 
estimates  of  the  deep  ocean  directional  spectrum. 

Work  during  the  past  year  was  concentrated  on  the  study  of  the 
characteristic  features  of  the  directional  spectra  at  Torrey  Pines  Beach. 
Particular  attention  was  paid  to  the  estimation  of  wave  energy  density  in  the 
shadowed  quadrants  of  the  local  directional  spectrum.  Three  directional  spectrum 
estimators,  two  of  which  were  developed  in  this  study,  were  compared  for  their 
quanitative  performance  in  the  definition  of  the  spectral  features.  The  results 
indicate  consistently  low  energy  density  levels  in  the  shadowed  quadrant 
associated  with  San  Clemente  Island  (SCI)  in  the  mid  frequency  range.  The 
spectra  at  higher  and  lower  wave  frequencies  occasionally  contain  significant 
energy  density  in  this  shadowed  quadrant. 

The  directional  sectors  a,,,  *  264.5-277°T  (SCI)  and  a0O>291.5°T  at 
Torrey  Pines  Beach  are  blocked  by  Islands  and  the  coastline.  The  island 
blocking/refraction  model  predicts  zero  energy  density  in  three  sectors  of 
the  local  deep  water  directional  spectrum.  Extensive  model  tests  were  used 
to  show  that  the  Maximum  Likelihood  Estimator  (MLE),  a  widely  used  directional 
spectrum  technique,  yields  biased-high  estimates  of  the  energy  density  in 


these  directional  "gaps".  Two  alternative  estimation  techniques  were 
developed.  The  first  method  attempts  to  deconvolve  the  smearing  effects 
due  to  the  limited  array  size  by  making  iterative  modifications  to  the  MLE 
estimate.  The  resultant  spectrum  (IMLE  estimate;  Pawka,  in  prep.)  smears 
through  the  array's  windowing  filter  to  yield  the  original  MLE  estimate.  In 
this  sense,  the  IMLE  estimate  is  a  possible  solution  for  the  "true"  direc¬ 
tional  spectrum.  The  second  estimator  attempts  to  "minimize  the  square 
window"  error  in  forming  its  estimate  at  each  look  angle.  The  derivation 
of  this  estimator  (MSWE)  follows  the  basic  design  of  the  Steered  Apriori 
Optimal  Estimator  developed  by  Regier  (1975).  The  directional  spectrum 
estimate  is  formed  from  weighted  sums  of  the  measured  cross-spectra.  The 
weights  are  adjusted  to  minimize  the  square  error  due  to  the  windowing  of 
the  measured  MLE  directional  spectrum  subject  to  the  constraint  that  energy 
density  at  the  look  angle  is  passed  unaltered.  A  noise  rejection  term  is 
also  added  to  the  minimization  equation.  A  discussion  of  the  development  of 
the  IMLE  and  MSWE  estimators  as  well  as  their  possible  applications  is  being 
prepared  for  publication. 

A  quantification  of  the  SCI  island  gap  definition  is  achieved  by 
normalizing  the  energy  density  at  the  center  of  the  gap,  Eg,  by  the  total 
frequency  band  energy,  Etp(f).  The  various  estimates  of  this  ratio  for  a 
data  set  on  5  March.. 1977  are  shown  in  Figure  Al-4.  The  MLE  estimates  of 
"normalized"  gap  energy  are  biased-high  relative  to  the  other  results  for 
all  wave  frequencies.  The  IMLE  and  MSWE  estimates  agree  well,  particularly 
in  the  range  0.121-0.145  Hz.  The  IMLE  estimates  in  the  mid  frequency  range 
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Figure  Al-4.  Normalized  San  Clemente  Island  gap  energy 
*-  «  density  on  5  March  1977.  The  results  are 

averaged  over  14.3  hours  of  continuous  data 
The  directional  bandwidth  is  1°  and  the  fre¬ 
quency  bandwidth  is  0.0078  Hz. 


are  more  stable  than  the  MSWE  results  but  show  a  small  relative  bias 
(high).  The  IMLE  directional  spectra  show  the  sharpest  definition  of  the 
gap  in  the  lowest  frequency  bands.  All  of  these  low  to  mid  frequency  results 
on  5  March  1977  data  match  predictions  made  from  model  testing  the  estimator 
responses  to  model  spectra  with  very  low  gap  energy  density. 

High  quality  data  from  10  days  of  the  Seasat  West  Coast  Experiment 
were  analyzed  with  the  MSWE  method.  The  results  in  the  frequency  range 
0.082-0.114  Hz  show  negligible  levels  of  energy  density  in  the  center  of  the 
SCI  gap.  The  average  directional  spectrum,  refracted  out  into  local  deep 
water,  for  0.114  Hz  (Figure  Al-5)  has  a  sharp  trough  at  ~271°T  which  corre¬ 
sponds  to  the  SCI  gap.  Variations  in  the  gap  energy  density  level  over  the 
data  runs  did  not  correlate  with  any  tested  parameter  of  the  wind  or  wave 
fields.  This  low  level  of  mean  gap  energy  density  is  consistent  with  the 
island  blocking/refraction  modeling  of  the  wave  transformations  from  deep 
ocean  to  the  coast. 

There  is  evidence  of  significant  SCI  gap  energy  levels  in  both 
the  low  and  high  frequency  regimes  for  this  data  set.  The  average  0.145  Hz 
spectrum,  also  shown  in  Figure  Al-5,  has  a  raised  gap  level  relative  to  the 
0.114  Hz  spectrum.  The  high  frequency  gap  energy  density  correlates  well 
with  the  local  wind  measured  at  the  NUC  tower,  Mission  Beach.  The  likely 
cause  of  the  low  frequency  gap  energy  density  was  not  determined.  The  maxi¬ 
mum  value  of  the  gap  energy  density  at  all  frequencies  was  201  of  the  peak 

% 

levels  in  the  directional  spectrum.  Therefore,  the  total  energy  contained 
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in  this  12.5°  sector  is  negligible  for  most  spectral  applications. 

There  is  a  15°  directional  sector  at  Torrey  Pines  Beach  centered 
on  305°T  which  is  totally  blocked  from  the  deep  ocean  but  has  a  fetch  of 
250-300  km.  The  MSWE  spectra  at  frequencies  greater  than  0.114  Hz  occa¬ 
sionally  show  a  directional  peak  associated  with  this  "inter-island"  fetch 
(Figure  Al-4).  The  energy  contained  in  this  spectral  peak  correlates  well 
with  the  local  wind.  The  peak  frequency  of  the  energy  associated  with  this 
directional  sector  during  local  wind  events  (wind  speeds  8-10  m/sec)  was 
0.137  Hz.  This  corresponds  well  with  the  values  computed  from  limited  fetch 
relationships  (Hasselmann  et  al . ,  1973).  This  peak  can  contain  up  to  30%  of 
the  total  energy  density  in  the  high  frequency  directional  spectra.  This 
significant  energy  level  combined  with  the  relatively  high  angle  of  these 
waves  with  respect  to  the  coast  cause  this  spectral  feature  to  be  an  important 
consideration  for  nearshore  processes  such  as  longshore  currents.  Therefore, 
local  wind  effects  have  to  be  taken  into  account  to  adequately  describe  the 
local  wave  field  for  wave  frequencies  greater  than  0.12  Hz. 

There  are  two  persistent  peaks  in  the  directional  spectra  at  Torrey 
Pines  Beach  associated  with  waves  propagating  to  the  site  from  the  North 
Pacific.  The  peak  with  mode  direction  '285-290°T  is  due  to  waves  coming 
through  the  primary  window  between  the  Channel  Islands  while  the  southern 
peak  (cu»~255-260°T)  is  due  to  refraction  over  banks  and  shoals  south  of  SCI 
(Pawka,  in  prep.).  Accurate  estimates  of  the  relative  energy  in  these  two 
prominant  peaks  is  a  necessity  for  quanitative  utility  of  the  directional 
spectra.  Ratios  of  the  energy  contained  in  these  two  peaks,  obtained  from 
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the  various  directional  spectrum  estimators  for  5  March  1977,  are  shown 
in  Figure  Al-6.  The  MSWE  and  IMLE  estimates  are  in  close  agreement  for 
the  frequency  range  0.114-0.145  Hz,  while  the  MLE  estimates  show  a  consider¬ 
able  bias.  The  MLE  estimates  have  the  most  significant  deviation  from  the 
results  of  the  other  methods  for  the  high  frequency  range. 

Data  comparisons  presented  here  and  extensive  model  testing  (Pawka, 
in  prep.)  have  shown  that  the  MLE  method  is  quantitatively  deficient  for  many 
desirable  applications  of  the  directional  spectra.  The  results  also  indicate 
improved  relative  performance  on  the  IMLE  and  MSWE  methods.  However,  no 
single  method  is  recommended  at  this  time  as  an  "optimal"  estimator  for 
general  application.  Intercomparison  of  the  results  of  several  methods  has 
proven  to  be  a  useful  analysis  procedure.  Additionally,  alternative  estima¬ 
tion  techniques,  tailored  to  address  specific  quanitative  questions,  may 
prove  to  be  more  desirable  in  limited  application. 
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Figure  Al-6.  Ratio  of  the  banks  quadrant  (250-269°T)  to  northern  window  (270-296°T) 
energy  for  the  various  directional  spectrum  estimators.  The  low 
frequency  ratio  Is  the  highest  due  to  the  strength  of  the  refraction 
process.  The  rise  of  the  relative  importance  of  the  banks  quadrant 
energy  at  higher  frequency  is  due  to  the  broadening  of  the  deep 
ocean  directional  spectrum. 


The  aim  of  the  shoaling  waves  task  has  been  to  develop  and  test 
a  model  describing  the  changes  undergone  by  wind  generated  surface  waves 
(4-18  sec  period)  as  a  broad  spectrum  of  such  waves  propogates  shoreward 
over  a  sloping  bottom.  A  one-dimensional  theoretical  model  based  on  the 
Boussinesq  equations  and  incorporating  the  physics  of  multiple  triad  near¬ 
resonance  has  been  derived  and  implemented  numerically.  The  model,  which 
assumes  that  all  waves  are  normally  incident  to  the  beach,  has  no  "free" 
or  empirically  determined  parameters.  A  single  field  experiment  can  thus 
be  used  to  determine  the  operational  validity  of  the  model  and  the  assump¬ 
tions  that  are  incorporated  in  its  derivations. 

Such  a  field  experiment,  involving  dense  instrumentation  of  the 
shoaling  region  from  10  m  depth  to  3  m  depth  was  conducted  from  June  to 
September  1980.  A  plan  view  of  the  experiment  is  shown  in  Figure  Al-7. 

Much  of  the  data  acquired  in  the  experiment  has  been  analyzed  in  light  of 
the  Boussinesq  model,  and  the  results  are  very  encouraging. 

As  shown  in  Figures  Al-8,  9,  10,  the  model  appears  to  do  an 
excellent  job  predicting  the  evolution  of  the  spectrum  of  sea  surface 
elevation  throughout  the  shoaling  region.  The  model  predicts  spectral 
shapes  significantly  better  than  does  linear  theory  over  a  broad  range  of 
input  conditions. 

Since  the  Boussinesq  model  predicts  component  amplitudes  and 
phases,  not  just  gross  spectral  quantities,  higher  order  moments  such  as 
skewness  of  depth  averaged  velocity  can  be  computed  and  compared  with  data. 
Both  the  model  and  the  data  show  a  significant  increase  in  negative  skewness 
(corresponding  to  higher  shoreward  velocities  and  peaked  wave  crests)  as  the 
waves  propogate  shoreward  through  the  shoaling  region,  as  shown  in 


21 


SHOALING  WAVES  EXPERIMENT 
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Figure  Al-7.  Plan  view  of  the  instrument  locations  in  the  1980 

shoaling  waves  field  experiment.  Approximate  depths 
are  given  near  the  right  side  of  the  figure. 
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Mgure  Al-8.  Comparison  plots  of  spectral  density  (cm/Hz)  vs.  frequency  (Hz)  at  6  on-offshore  locations. 
Each  plot  shows  the  actual  data  at  the  Indicated  depth,  the  prediction  of  the  Bousslnesq  model,  and  the 
prediction  of  linear  theory.  Both  the  Boi'v*; inesq  and  linear  theory  predictions  are  based  on  a  240  mode 
decomposition  of  the  10  m  data.  This  day  ?*.  Sept  80)  has  very  low  energy  and  a  broad  spectrum.  Note  that 
both  linear  theory  and  the  Boussinesq  mod,.;  predict  the  spectrum  well. 
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ame  as  Figure  Al-8,  but  this  day  (11  Sept  80)  Is  characterized  by  a  very  narrow 
pectrum  as  measured  in  10  m  depth.  Linear  theory  cannot  predict  the  increase 
f  harmonic  energy  as  the  waves  propogate  shoreward. 
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Figure  Al-11.  Although  quantitative  agreement  between  model  and  data  is 
not  as  good  as  for  spectral  comparisons,  the  model  predictions  are  always 
within  a  factor  of  2  of  the  data. 

Dissipation  of  wave  energy  through  the  shoaling  region  is  expected 
if  energy  is  being  extracted  from  the  waves  due  to  bottom  stress.  Both  the 
data  and  the  Boussinesq  model  exhibit  non-conservation  of  energy  flux  com¬ 
puted  from  linear  theory  through  the  shoaling  region  (Figure  Al-12).  As 
the  model  contains  neither  dissipation  nor  generation,  the  apparent  decreases 
and  increases  of  linear  energy  flux  are  an  indication  of  the  inappropriate¬ 
ness  of  linear  ideas  applied  to  even  weakly  nonlinear  systems.  The  fact 
that  the  model  predictions  are  quantitatively  similar  to  data  at  the  3% 
level  indicate  that  dissipation  in  the  shoaling  region  before  the  break 
point  is,  in  fact,  small. 

In  a  model  independent  analysis,  the  data  obtained  by  the  two  wave 
directional  arrays  (enclosed  by  dashed  lines  in  Figure  Al-7)  is  being  used 
to  determine  the  amount  and  frequency  distribution  of  seaward  propogating 
(reflected)  energy  in  the  wind  wave  band.  Preliminary  results  show  that 
the  outgoing  energy  is  small  (<10%)  compared  to  the  incoming  energy  over  a 
wide  range  of  conditions. 

Theoretical  model.  In  order  to  be  applicable  to  field  situations, 
any  theoretical  shoaling  model  must  allow  for  a  complicated  wave  field  com¬ 
posed  of  a  braod,  arbitrary  frequency  spectrum. 

Two  basic  types  of  models  have  been  used  in  the  past  to  predict 
the  wave  field  parameters  of  surface  elevation,  pressure,  and  particle  velo¬ 
city  in  the  shoaling  region:  Linear  theory,  and  forced  nonlinear  (Stokes) 
theories.  Linear  theory  is  based  on  the  assumption  that  nonlinear  terms  in 
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Figure  Al-12.  Comparison  plot  of  R(h,hQ)  versus  depth  for  11  September  80  (hQ=10.3m).  If 
linear  energy  flux  were  conserved,  R(h,hQ)  would  be  a  constant  R ( h , hQ }  =  1.  The  fact  that 
the  model  compares  well  with  data  indicates  that  R(h,ho)  is  not  a  conserved  quantity  for  a 
nonlinear  system  such  as  shoaling  surface  gravity  waves. 


inviscid,  irrotational  equations  of  motion  and  boundary  conditions  are 
unimportant.  The  resulting  equations  can  be  solved  analytically  even 
for  a  wave  field  of  arbitrary  spectral  shape.  The  solutions  conserve 
energy  flux  and  predict  a  general  increase  in  sea  surface  elevation  as 
the  waves  propogate  shoreward.  However,  the  model  has  many  drawbacks: 
bottom  slopes  must  be  extremely  small,  and  amplitudes  of  waves  long  in 
comparison  to  the  depth  must  be  unrealistically  tiny  for  the  linearization 
to  remain  consistent.  Since  the  theory  itself  predicts  increasing  ampli¬ 
tude  for  decreasing  depth,  this  latter  condition,  even  if  true  initially, 
is  quickly  violated.  In  addition,  linear  theory  predicts  that  wave  forms 
remain  symmetric  throughout  the  shoaling  region  and  cross-spectral  energy 
transfers  do  not  take  place,  both  in  obvious  contradition  to  observations. 

Forced  nonlinear  theories  take  into  account  the  modifications  to 
a  single  linear  wave  train  when  its  amplitude  becomes  finite.  The  solution 
for  the  wave  train  at  any  given  order  is  composed  of  the  fundamental  wave 
and  its  forced  harmonics.  Considerable  effort  has  been  expended  in  the 
development  and  testing  of  such  theories.  Well  known  perturbation  methods 
can  be  used  to  derive  algebraic  expressions  for  the  relative  amplitudes  of 
the  fundamental  and  harmonics.  Flick,  et  al  (1981)  have  shown  that  such 
solutions  valid  in  deep  water  can  be  matched  to  cnoidal  wave  solutions 
valid  in  shallow  water.  In  addition,  these  workers  under  ONR  sponsorship 
have  conducted  extensive  laboratory  testing  of  the  solutions.  Unfortunately, 
monochromatic  forced  solutions  are  inappropriate  to  cases  where  the  offshore 
wave  spectrum  is  broadbanded  in  frequency.  Algebraic  complexity  grows 
rapidly  as  more  fundamental  waves  are  added  and  the  matching  to  shallow 
water  solutions  is  no  longer  valid.  Even  in  the  monochromatic  case,  the 


model  cannot  predict  the  evolution  of  the  assymmetric  wave  forms  charac¬ 
teristic  of  waves  near  breaking. 

In  response  to  the  above  deficiencies  of  both  linear  theory  and 
forced  nonlinear  theories,  we  have  developed  a  shoaling  model  based  on  the 
Boussinesq  equations  and  the  physics  of  multiple  triad  near  resonance  across 
the  entire  wind  wave  frequency  band.  The  "Boussinesq  model"  currently  allows 
for  waves  propogating  only  normal  to  a  beach  with  parallel  contours.  In 
addition,  it  is  not  clear  how  the  model  can  be  matched  in  an  analytic  way 
to  Stokes-type  descriptions  of  deep  water  waves.  However,  the  model  pre¬ 
dicts  the  evolution  of  both  broad  and  narrow  banded  input  spectra,  and  it 
allows  for  the  development  of  asymmetric  wave  forms.  For  a  detailed  deri¬ 
vation  of  the  one  dimensional  Boussinesq  shoaling  model,  the  reader  is 
referred  to  Task  (80-2)  of  the  1980  ONR  progress  report. 

The  model  is  cast  as  a  set  of  coupled,  first  order  differential 
equations  for  the  on-offshore  evolution  of  amplitudes  and  phases  of  the 
Fourier  modes  of  the  wave  field.  The  model  contains  no  adjustable  parameters 
or  empirically  determined  constants.  It  can  thus  be  directly  tested  against 
data.  The  plan  view  of  an  experiment  designed  to  test  the  model  is  shown 
in  Figure  Al-7.  The  experiment  design  was  discussed  extensively  in 
Task  (81-2)  of  the  1981  ONR  proposal.  Comparison  between  data  from  the 
on-offshore  instrument  transect  and  predictions  of  the  Boussinesq  model  is 
straight  forward.  Each  of  the  17.1  minute  records  of  data  at  each  on-offshore 
location  are  reduced  and  transformed  to  yield  Fourier  components  of  depth 
averaged  velocity.  A  spectrum  composed  of  240  modes  from  location  P6 
(approximately  10  m  depth)  Is  used  as  the  Initial  condition  for  the  numeri¬ 
cal  integration  of  the  model's  coupled  evolution  equations. 
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Spectral  Evolution.  Spectral  comparisons  between  the  model  pre¬ 
dictions  and  data  can  be  made  at  each  on-offshore  location.  Figure  Al-8 
shows  such  comparisons  at  6  locations,  ranging  in  depth  from  9. 0-3. 5  m 
on  5  Sept  80.  The  dashed  third  line  on  each  plot  represents  linear  theory 
predictions  of  spectral  density.  This  particular  data  is  an  example  of  a 
broad  banded,  very  low  energy  wave  field  (Hs  in  9.9  m  depth  is  38  cm). 

As  expected  with  such  small  waves,  nonlinear  effects  are  small,  and  both 
linear  theory  and  the  Boussinesq  model  predict  spectral  shape  well. 

Figure  Al-9  is  a  similar  comparison  on  11  Sept  80.  The  offshore 
spectrum  is  quite  narrow,  with  most  wave  energy  concentrated  in  a  peak  cen¬ 
tered  at  0.065  Hz.  (Hs  is  90  cm).  Spectral  evolution  is  seen  to  be 
characterized  by  the  emergence  of  a  major  narrow  peak  centered  at  0.127  Hz, 
and  a  smaller  peak  at  0.190  Hz.  Although  linear  theory  does  not  predict 
these  peaks  at  all,  the  Boussinesq  model  accurately  predicts  both  their 
locations  and  amplitudes  throughout  the  shoaling  region. 

Figure  Al-10  presents  yet  a  third  comparison,  using  data  obtained 
on  9  Sept  80.  In  this  case,  the  spectrum  in  10.3  m  depth  is  fairly  broad, 
and  the  waves  are  of  moderate  height  (H$  is  66  cm).  A  dramatic  evolution 
of  the  spectrum  is  observed,  with  the  entire  high  frequency  region  increas¬ 
ing  in  energy.  Once  again,  linear  theory  does  not  properly  predict  the 
spectral  evolution,  while  the  Boussinesq  model  compares  very  well  with  the 
data.  The  Boussinesq  model  is  thus  seen  to  predict  spectral  evolution 
well  over  a  wide  variety  of  wave  conditions  and  spectral  shapes. 

Skewness.  Third  and  higher  moments  of  the  wave  field  have  been 
shown  to  be  important  to  sediment  transport  processes  (Task  B2).  In  water 
deeper  than  10  m,  the  wave  field  is  to  a  good  approximation  Gaussian,  and 
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thus  all  odd  moments  are  exceedingly  small.  Previous  laboratory  and 
scattered  field  measurements  indicated  that  skewness  of  surface  elevation 
or  depth  averaged  velocity  increases  in  a  negative  way  (corresponding  to 
sharp  wave  crests  and  large  maximum  onshore  fluid  velocities).  To  our 
knowledge,  this  evolution  has  never  been  measured  in  a  detailed  way  in 
the  field.  The  Fourier  coefficients  of  the  data  at  each  on-offshore 
location  can  be  used  to  calculate  velocity,  and  the  results  for  9  Sept  80 
are  presented  in  figure  Al-11.  Also  shown  is  the  skewness  as  calculated 
from  the  Boussinesq  model.  Although  the  model  skewnesses  are  everywhere 
larger  (in  a  negative  sense)  than  those  calculated  from  data,  the  model 
predictions  are  well  within  a  factor  of  2  of  the  data.  This  appears  to 
be  a  general  result  for  all  data  analyzed  so  far. 

Dissipation.  The  amount  of  energy  lost  from  the  wave  field  is 
of  crucial  importance  to  studies  of  sediment  transport  outside  the  break 
point.  Since  wave  breaking  and  associated  turbulent  dissipation  is  not  a 
major  factor  in  the  shoaling  region  (in  comparison  to  the  surf  zone), 
detailed  measurements  of  the  "dissipation"  of  wave  energy  in  this  region 
could  provide  a  realistic  upper  bound  for  the  energy  available  to  move  sedi¬ 
ments.  The  Boussinesq  shoaling  model  alone  provides  no  information  on  dis¬ 
sipation,  as  it  was  assumed  at  the  outset  that  all  viscous  terms  in  the  full 
Navier-Stokes  equations  could  be  neglected. 

The  model  can  be  used  however,  to  show  that  energy  flux  as  calcu¬ 
lated  by  linear  theory  is  not  conserved  over  distances  comparable  to  the 
on-offshore  extent  of  the  shoaling  region.  The  mathematics  is  lengthy  and 
the  formal  ordering  arguments  needed  In  even  the  simplest  cases  are  subtle, 
and  will  not  be  reproduced  here.  We  present  Instead  comparisons  between 
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total  linear  energy  flux  as  calculated  from  the  data  at  various  on-offshore 
locations,  and  the  same  quantity  as  calculated  from  model  integrations.  At 
any  depth  or  on-offshore  location,  the  total  linear  energy  flux  density  can 
be  defined  as: 

F(h)  =/[a(o)]2  ||  do 

where  [a (o) ]  is  the  spectral  density  of  surface  elevation  at  frequency 
o,  and  is  "group  velocity"  calculated  from  the  linear  dispersion  rela¬ 
tion.  In  a  dissipationless  linear  system,  F  is  constant,  independent  of 
h  or  on-offshore  distance.  A  measure  of  non-conservation  of  F  is  given 
by: 

R(h,h0)  =  F(h)/F(h0) 

The  effect  of  dissipation  on  the  linear  system  would  be  to  cause  R(h,hQ) 
to  decrease  monotonically  for  h<hQ  for  waves  propogating  towards  the 
shore.  An  observed  increase  in  R(h,hQ)  could  be  due  to  a  localized  source 
of  wave  energy,  which  is  unlikely. 

The  solid  line  in  Figure  Al-12  presents  a  typical  example  of  the 
behavior  of  R(h,h0)  as  calculated  from  the  data.  In  this  case,  hQis  10.3  m. 
Note  that  R(h,hQ)  both  increases  and  decreases  appreciably  through  the 
shoaling  region.  The  dashed  line  presents  an  identical  calculation  made 
using  the  results  of  a  model  integration.  Here  also  R(h,hQ)  clearly  is 
neither  constant  nor  is  its  behavior  monotonic.  It  should  be  remembered 
that  the  model  lacks  any  mechanism  for  either  dissipation  or  generation. 
Thus,  the  nonconservation  of  F(h)  in  this  case  is  a  strong  indication  that 
F(h)  is  not  a  conserved  quantity  for  a  weakly  nonlinear,  multiply  resonant 
system  such  as  shoaling  waves.  The  fact  that  model  predictions  of  linear 


energy  flux  are  quantitatively  similar  (±3%)  to  local  calculations  from 
the  data,  can  be  used  to  construct  the  argument  that  total  dissipation 
from  the  wind  wave  frequency  band  must  be  small  through  this  shoaling 
region.  It  is  clear,  because  of  nonlinear  effects,  that  experiments 
which  measure  only  wave  field  parameters  have  limited  utility  as  direct 
tests  of  the  importance  of  dissipation. 

Outgoing  Energy.  Data  from  the  two  directional  arrays  (enclosed 
by  dashed  lines  in  Figure  Al-7  )  can  be  used  to  construct  full,  360°  direc¬ 
tional  spectra  in  the  wind  wave  frequency  band.  It  is  hoped  that  this  data 
will  eliminate  the  need  for  doing  a  completely  new  experiment  to  test  a 
future  two-dimensional  shoaling  model.  In  a  purely  descriptive  sense, 
the  data  from  the  arrays  can  be  used  to  characterize  outgoing  (seaward 
propogating)  energy  in  the  shoaling  region. 

Such  energy  could  be  caused  either  by  gradual  reflection  from  the 
sloping  bottom,  or  the  highly  nonlinear  processes  dominating  the  surf  zone. 
Preliminary  results  from  the  current  experiment  indicate  that  outgoing 
energy  is  small  (  10%  of  incoming  energy)  in  the  wind  wave  band.  The 
Maximum  Likelihood  Estimator  used  to  process  array  data  cuases  energy  to 
spuriously  appear  to  be  seaward  propogating.  Thus  extensive  model  testing 
is  necessary  to  quantitatively  discuss  outgoing  energy.  However,  these 
preliminary  results  confirm  the  utility  of  linear  longshore  arrays  for 
measurements  of  nearshore  wave  directions. 


A2.  Shelf  Waves  (R.  T.  Guza) 

Shelf  waves  are  defined  here  as  oscillations  confined  to  coastal 
regions  because  of  topographic  trapping  effects  associated  with  shelf 
and  beach  bathymetry.  Shelf  waves  considered  in  this  task  have  period' 
ranging  from  1  to  10 1 s  of  minutes,  and  horizontal  length  scales  from 
100’s  of  meters  to  10's  of  kilometers.  Various  kinds  of  shelf  waves  are 
thought  to  be  important  elements  in  harbor  seiching,  currents  on  the 
continental  shelf  and  surf  beat.  The  long  term  goal  of  this  task  is  to 
assess  the  relative  importance  of  shelf  waves  in  various  physical  environ¬ 
ments  and  to  identify  the  mechanisms  leading  to  their  excitation.  Effort 
so  far  have  been  concentrated  on  surf  beat,  since  this  is  a  phenomenon  of 
first  order  importance  to  the  inner  surf  zone  velocity  field. 

Surf  beat  are  the  low  frequency  (periods  of  one  to  several  minutes) 
oscillations  which  typically  contain  a  substantial  portion  of  the  energy 
observed  in  inner  surf  zone  elevation  and  velocity  fields.  Surf  beat  has 
its  maximum  energy  levels  in  the  run-up.  During  the  past  contract  year, 
a  large  amount  of  run-up  data  consisting  of  time  histories  of  the 
swash/dry  beach  intersection  point  previously  collected  at  Torrey  Pines 
Beach  has  been  analyzed.  The  goal  of  the  analysis  is  to  empirically 
relate  the  level  of  surf  beat  activity  in  the  surf  zone  to  the  amount  of 
incident  wind-generated  wave  energy  measured  with  an  offshore  sensor. 

An  important  conclusion  is  that  the  magnitude  of  typical  swash 
excursions,  increases  approximately  linearly  with  increasing  incident 
wave  height  H$  .  Figure  A2-1  shows  the  linear  trend  (solid  line)  for  a 
variety  of  incident  wave  conditions.  When  the  swash  excursions  across  the 


beach  face  are  transformed  to  vertical  displacements,  significant  vertical 
run-up  heights  are  found  to  be  roughly  70%  of  the  significant  height  of 
the  incident  wind  waves  at  10  meter  depth. 

Previous  ONR  supported  work  (for  example  Inman  et  al . ,  1968)  pointed 
out  that  typical  run-up  spectra  show  low  frequency  (surf  beat)  oscillations 
to  be  more  energetic  than  run-up  motions  at  wind  wave  frequencies.  The 
present  data  set  was  used  to  quantify  that  observation  by  dividing  the 
run-up  energy  into  low  and  high  frequency  components.  A  frequency  of 
0.05  Hz  (period  20  s)  was  chosen  as  the  highest  frequency  considered  to 
be  part  of  the  surf  beat  motions.  Figure  A2-2  shows  that  the  high  fre¬ 
quency  portion  of  the  run-up  (numbers  in  circles)  remain  approximately 
constant  with  increasing  Hs>  while  the  surf  beat  varies  linearly.  On 
the  day  with  the  largest  incident  waves,  surf  beat  swash  amplitudes  are 
about  six  times  that  of  the  wind  wave  frequencies.  In  the  inner  surf 
zone,  the  presence  of  large  amplitude  surf  beat  motions  is  probably  the 
most  important  dynamic  consequence  of  energetic  incident  waves  and  a 
wide  surf  zone.  Details  of  the  above  work  are  given  in  Guza  and 
Thornton  (in  press). 

The  motions  at  surf  beat  periods  are  probably  related  to  variations 

in  the  heights  of  incident  waves.  There  is  also  a  surf  zone  response  to 

the  mean  incident  wave  height.  The  mean  water  line  at  the  shoreline 

is  elevated  ("set-up")  above  its  location  in  the  absence  of  waves.  The 

Torrey  Pines  run-up  data  (Figure  A2~3)  shows  that  the  shoreline  set-up 

n  is  about  17%  of  the  deep  water  significant  wind  wave  height,  HSoo 
M 

(Guza  and  Thornton,  1981). 


Figure  A2-3.  Shoreline  set- 
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A3.  Surf  Zone  Dynamics 

3.1  Turbulence  (R.  E.  Flick  and  R.  T.  Guza) 

Dividing  surf  zone  fluid  dynamics  in  a  frequency/wavelength 
sense,  a  substantial  amount  is  known  about  motions  in  the  fluid  interior 
in  the  typical  wind  wave  frequency  band,  less  about  surf  beat,  still  less 
about  the  high  frequency  gravity  wave  range  and,  essentially  nothing  about 
turbulent  fluctuations.  In  spatial  terms,  the  least  understood  areas  are 
very  near  the  water  surface  and  bottom.  Unfortunately,  the  single  most 
important  fluid  dynamical  topic  for  sediment  transport  is  probably  the 
turbulent  bottom  boundary  layer,  as  detailed  in  Section  Bl. 

Work  on  this  task  has  begun  during  the  past  contract  year  and 
it  is  closely  related  to  the  bedload  investigation  of  oscillatory  bursting 
discussed  in  Section  Bl.l.  This  task  proposed  a  field  study  of  the  gross 
characteristics  of  naturally  occurring,  rough  turbulent  flow  in  and  near 
the  surf  zone  under  sheet  flow  conditions. 

Progress  has  been  made  in  the  design  and  fabrication  of  temperature 
and  velocity  sensing  instruments  and  mounting  platform's  for  their  deployment. 
A  rigid  platform  has  evolved  based  on  the  very  strong  and  extensively 
tested  "surf  spider".  This  triangle  based,  slender  legged  apparatus  was 
developed  to  hold  current  meters,  wave  staffs  and  other  instruments  in  the 
swash  and  shallow  surf  zone.  The  legs  were  modified  to  hold  the  spider 
rigid  underwater  in  such  a  way  to  result  in  minimum  disturbance  to  the  flow. 

"fc 

A  cone  probe  hot  film  anemometer  sensor  has  been  purchased  and 
packaged  for  deployment.  Laboratory  calibrations  for  speed  and  directional 
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sensitivity  are  underway.  The  low-frequency  portion  of  the  anemometer 
signal  will  be  compared  to  an  electromagnetic  current  meter  during 
actual  data  runs  to  assess  calibration  reliability  of  the  cone  probe. 

Fast  response  thermistor  probes  and  circuitry  have  also  been 
assembled  and  are  awaiting  calibration.  This  system  will  be  used  to 
detect  small  scale  temperature  fluctuations  which  may  affect  the  hot- 
film  anemometer  measurements. 
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A4.  Current  Measurements  by  Geomagnetic  Induction  (R.  E.  Flick) 

Progress  during  the  contract  year  has  been  made  in  the  design 
and  construction  of  stable  silver-silver  chloride  electrodes  and  of  a 
salt  bridge  water-chopped  system.  Both  these  will  be  deployed  in  early 
fall  1981  across  Mission  Bay  Channel.  The  goal  of  this  deployment,  de¬ 
tailed  in  last  year's  proposal,  is  to  measure  the  tidal  transport 
induced  electric  potential  across  the  channel.  Corrections  for  bottom 
conductivity  will  then  be  applied  and  the  resulting  estimates  of  volume 
transport  compared  with  predictions  from  various  theories  which  seek  to 
model  the  flow  forcing.  To  the  extent  that  long  waves  on  the  shelf  force 
bay  mode  resonances,  this  task  is  connected  to  exploratory  measurements 
underway  and  proposed  under  Task  A2. 

The  electrodes  are  much  larger  than  those  previously  used  in  the 
unsuccessful  attempts  to  measure  cross-shore  electric  potentials  (Flick  and 
Parker,  1980).  They  are  also  more  stable  and  less  sensitive  to  pressure 
fluctuations.  These  "quarter-pounder"  electrodes  have  been  developed  with 
the  advice  and  cooperation  of  Dr.  Jean  Filloux  of  Scripps. 

The  salt  bridge  water-chopped  instrumentation  has  been  designed 
to  act  like  a  double  pole-double  throw  switch  which  alternately  connects 
each  of  the  two  electrodes  to  the  opposite  side  of  the  flow.  This  must  be 
done  to  avoid  polarization  and  compensate  for  drift  of  the  electrodes  over 
the  planned  one  to  two  month  measurement  time. 


B.  FLUID-SEDIMENT  INTERACTIONS 
"BED  RESPONSE" 

The  ultimate  objective  of  this  contract  is  the  prediction  of 
shelf  and  beach  forms  given  the  driving  forces,  the  general  shelf 
bathymetry  and  the  nature  of  the  bottom  sediments.  Accordingly,  the 
tasks  under  A  are  directed  towards  defining  the  relevant  physical 
properties  of  the  driving  forces;  while  the  tasks  under  B,  the  bed 
response  to  fluid-sediment  interactions,  lead  to  an  understanding 
of  the  shelf  and  beach  forms  that  are  the  central  theme  of  this  con¬ 
tract.  The  tasks  under  C  are  largely  exploratory  in  nature,  with  the 
general  objective  of  determining  the  extent  to  which  signatures  from 
remote  sensors  can  be  used  to  define  the  broader  spacial  character¬ 
istics  of  driving  forces,  the  larger  scale  ocean  surface  patterns 
that  are  characteristic  of  coastal  zone  dynamics,  and  the  large  scale, 
longer  term  changes  that  occur  in  shoreline  configuration. 

All  bed  forms  result  from  fluid-sediment  interactions  that 
genetically  divide  into  two  primary  bed  form  classes,  "forced  forms", 
and  "interactive  forms";  or,  a  combination  of  these  two  (Task  B3.1). 

Since  the  generality  of  this  bed  ‘orm  classification  appears 
to  be  valid,  it  forms  a  valuable  working  hypothesis  for  the 
study  of  fluid-sediment  interactions.  However,  we  know  little 
about  the  mechanics  of  the  interactions,  nor  the  nature  and  strength 
of  the  feedback  mechanisms.  Since  the  bed  responds  to  the  ^luid  forces 
by  changing  its  form  and  the  bed  form  in  turn  feeds  momentum  and  vor- 
ticity  back  into  the  fluid,  interactive  phenomena  evolve  that  are  far 
more  complex  than  the  relatively  simple  problem  of  a  fluid  stress  over 
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a  rough  but  immovable  bed.  In  fact,  the  granular  bed  itself  responds 
to  stress  by  expanding  its  volume  (dil itation)  so  that  the  nature  of  the 
granular  motions  within  the  bed  become  functions  of  both  the  external 
fluid  stress  and  the  volume  concentration  of  the  solid  particles  in  the 
bed.  Althoughthis  property  of  granular-fluid  media  is  fundamental  to 
bed  load  transport  and  to  bed  forms,  its  mechanics  remain  virtually 
unknown  (Bagnold,  1954;  Inman,  1963;  Bailard  and  Inman,  in  prep; 

Savage,  1979). 

As  a  consequence,  three  tasks,  Bl.l  through  B1.3,  under  this 
section  are  studies  of  sand  bed  response  to  fluid  stress.  Task  Bl.l  is 
a  study  of  bed  response  to  near-breaking  waves  and  an  investigation  of 
the  possibility  that  this  intense  shear  involves  "oscillatory"  bursting. 
Task  B1.2  is  a  study  of  the  rheology  of  granular-fluid  material,  a  neces¬ 
sary  prerequisite  to  understanding  the  quantitative  aspects  of  sediment 
transport.  Task  B1.3  is  to  deploy  a  bedload  meter  to  measure  the  thick¬ 
ness  and  velocity  of  the  bedload  under  waves  and  currents  in  the  field. 
Once  the  bedload  due  to  waves  and  currents  is  known  at  least  empirically, 
then  a  rational  approach  to  on-offshore  transport  and  the  prediction  of 
beach  slopes  can  be  employed  (Task  B2). 

Task  B3  is  a  fundamental  study  of  the  formation  of  shelf  and 
beach  features  which  will  supplement  the  study  of  "Bed  forms  produced  by 
fluid-granular  interaction"  begun  as  ONR  Task  (76-4).  This  general  task 
includes  the  fundamental  study  "Bedforms  produced  by  fluid-granular 
interaction"  (B3.1),  and  a  definitive  study  of  the  limiting  amplitude  of 
beach  cusps  (B3.2).  Task  B4  will  provide  information  on  coastal  forms 
and  the  budget  of  sediment  in  tideless  seas  and  may  be  expanded  to  study 
the  mechanism  of  wave  attenuation  in  shallow  water. 


Bl.l  "Oscillatory  Bursting"  (D.  L.  Inman  and  S.  A.  Jenkins) 

Bursting  is  a  high  Reynolds  number  mechanism  in  both  steady  and 
unsteady  shear  flow  for  advecting  vorticity  from  the  bottom  boundary 
layer  into  the  interior  of  the  fluid.  Vorticity  leads  to  turbulent 
vertical  velocities  which  in  turn  suspends  sediments.  The  burst  cycle 
starts  with  an  instability  in  the  boundary  layer  which  subsequently 
degenerates  under  certain  conditions.  The  instability  could  result 
from  an  inflection  in  the  velocity  profile  near  the  boundary  due  to 
adverse  pressure  gradients  from  certain  types  of  bed  forms,  or  from  a 
flexible  or  permeable  boundary,  or  even  from  external  pressures  such  as 
gravity  waves.  Vorticity  is  concentrated  in  the  lee  of  the  instability, 
and  accumulates  as  small  "turbulent  spots"  (e.g.  Schlichting,  1979,  Ch  17). 
Inertial  shear  waves  known  as  "Tollmien-Schlichting  (T-S)  waves  form  in  the 
shear  flow  in  the  upper  portion  of  the  boundary  layer  over  the  curvature 
formed  by  the  turbulence  spots.  These  T-S  wav,es  grow  and  steepen  when 
viscous  stresses  are  no  longer  sufficient  to  damp  them  out.  This  insta¬ 
bility  occurs  when  the  wave  length  of  the  T-S  waves  become  smaller  than 
approximately  66  where  6  is  the  boundary  layer  thickness.  Upon  becom¬ 
ing  unstable  the  T-S  wave  breaks  or  bursts  into  the  interior  of  the  fluid. 
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Oscillatory  Bursting:  An  Hypothesis.  It  has  been  thought 
that  bursting  does  not  occur  under  oscillatory  waves.  Jackson  (1976, 
p.  537)  states:  "The  writer  knows  of  no  measurements  of  the  burst 
cycle  in  the  thin  boundary  layers  beneath  wind-generated  shoaling  waves. 

The  swiftly  reversing  flow  virtually  precludes  rational  measurement"  ... 
and  the  formation  of  burst  cycles.  Also,  inspection  of  the  character¬ 
istics  favoring  bursting,  suggest  that  under  oscillatory  motion  there  is 
much  less  likelihood  of  favorable  characteristics,  except  under  very 
special  circumstances  (Hinze,  1975). 

On  the  other  hand,  there  are  a  number  of  observations  that 
suggest  that  a  form  of  bursting  does  occur  under  shoaling  waves.  For 
example,  Inman  (1957,  p.  23)  observed  that  the  higher  velocities  associ¬ 
ated  with  waves  in  shallow  water  "appear  to  place  the  material  in  a  state 
of  'sub-suspension'  with  the  passage  of  each  wave  crest.  In  this  case 
the  bottom  material  is  moved  on  and  off  shore  as  a  thin  but  dense  layer 
of  suspended  particles,  and  the  discrete  cloud  of  suspended  particles 
which  form  at  the  crest  of  the  trochoidal  shaped  ripples  appears  to  be 
absent.  The  thickness  of  the  layer  of  'sub-suspension'  although  somewhat 
variable,  is  usually  of  the  order  of  a  few  centimeters,  and  appears  to 
depend  on  the  size  of  the  sand  and  on  the  orbital  velocity  and  periodi¬ 
city  of  the  waves." 


In  general  if  one  observes  sand  motion  through  a  face  plate 
during  the  passage  of  swell  waves  near  the  surf  zone,  the  following  sequence 
will  occur.  Assume  the  waves  are  of  10  to  15  second  period,  1  meter  height, 
over  a  fine  sand  bottom  and  travel  in  water  2  to  3  meters  deep 
(Figure  Bl-1).  As  the  wave  crest  approaches,  the  onshore  velocity  will 


Figure  Bl-1.  Schematic  phase  diagram  depicting  development  of  unstable  oscillatory 
boundary  layer  under  a  progressive  wave  In  shallow  water.  Note  that 
both  pressure  gradient  (dp/dx)  and  boundary  layer  velocity  (u£)  are 
unstable  following  crest  passage.  Also,  tne  thickness  of  the  unstable 
portion  of  boundary  layer  Increases  from  phase  cycle  0  through  2. 


rapidly  increase  to  a  threshold  value  and  sand  grains  will  begin  to 
move  as  bedload.  Then  suddenly,  just  as  the  crest  passes,  the  sand 
appears  to  "boil"  and  a  dense  cloud  of  sand  approximately  3  to  5  cm  thick 
moves  onshore  until  the  near  bed  motion  is  outrun  by  the  faster  traveling 
wave  crest.  The  sand  cloud  then  drops  back  to  the  bed  and  the  bed  is 
observed  to  be  nearly  flat,  with  just  a  faint  trace  of  residual  ripples. 
The  ripples  have  wavelengths  of  about  6-7  cm  and  heights  of  1/2  cm  or  less 
(e.g.  Dingier  and  Inman  1976,  fig.  9).  Upon  closer  examination  it  will  be 
noted  that  the  sand  "boils"  when  they  occur,  have  an  on-offshore  separa¬ 
tion  approximately  equal  to  the  residual  ripple  wavelength,  and  their 
densest  portions  have  more  along-crest  coherence  as  though  they  were  asso¬ 
ciated  with  the  remnant  ripple  crests. 

The  above  phenomena  is  thoughtbyus  to  represent  "oscillatory 
bursting"  associated  with  "vortex  spots"  in  the  oscillatory  boundary. 

It  is  thought  that  the  T-S  waves  burst  because,  unlike  the  potential 
oscillation  in  the  water  above,  the  vertically  dependent  phase  shift  of 
the  boundary  layer  oscillation  allows  it  to  decelerate  at  the  time  the 
pressure  gradient  turns  negative  just  as  the  wave  crest  passes 
(Figure  Bl-1).  Oscillatory  bursting  is  limited  to  a  few  centimeters 
from  the  bottom  because  only  a  limited  region  of  the  boundary  layer 
supports  decelerating  flow  under  negative  pressure  gradients  imposed 
by  the  wave  profile,  and  exists  only  for  a  fraction  of  a  wave  cycle. 
Therefore  the  penetration  distance  of  oscillatory  bursting  is  limited 
by  the  onset  of  flow  reversal. 
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The  laminar  boundary  layer  under  oscillatory  waves  is  given 


and  would  be  about  2  mm  and  3  mm  thick  for  waves  of  10  and  15  second  period 
respectively.  However,  Sleath  (1970)  shows  that  the  boundary  layer  over 
a  flat  sand  bottom  is  rough  and  can  be  represented  approximately  as 

6r  -  C  6 

where  c  is  a  dimensionless  number  that  is  a  function  of  the  Reynolds 
number  for  the  sand  grains  and  wave  boundary  layer.  Re  *  (umD/v) (D/6) . 

Over  sand  beds  c  is  of  order  10,  and  thus  the  actual  boundary  layer 
thickness  is  about  2  cm  and  3  cm  respectively  for  10  and  15  second  waves. 

This  boundary  layer  thickness  is  in  approximate  agreement  with  the 
observed  3-5  cm  height  of  the  sand  "boils"  associated  with  shallow 
water  waves. 

Experimental  Investigation.  The  above  hypothesis  for  "oscil¬ 
latory  bursting"  is  now  being  tested  by  field  observations  just  outside 
the  surf  zone,  using  relatively  simple  sensors,  and  a  short  run  with  a 
hot  film  anemometer  (a  relatively  complex,  fragile  sensor).  The  hot 
film  anemometer  will  be  used  in  cooperation  with  Task  A3.1. 

The  study  uses  a  motion  picture  camera,  a  pressure  sensor  feeding 
synchronous  timing  lights,  and  a  high  resolution  sonar  scanner  (Dingier  and 
Inman,  1976).  If  sand  "boiling"  is  found  to  coincide  with  nearcrest  passage 
and  to  have  the  proper  length  scales  as  postulated,  then  hot  film  anemometers 
will  be  used  to  sense  the  presence  of  T-S  waves. 


B1.2  Rheology  of  Granular-Fluid  Material  (D.  M.  Hanes  and  D.  L.  Inman) 


One  of  the  least  understood  parameters  influencing  bedload 
sediment  transport  is  the  magnitude  of  the  bedload  itself!  This  is  a 
manifestation  of  the  general  lack  of  understanding  of  the  physics 
involved  in  the  deformation  and  flow  of  granular  and  granular-fluid 
media. 

In  1954,  R.  A.  Bagnold  pioneered  the  field  of  granular-fluid 
mechanics  by  studying  the  flow  of  neutrally  bouyant  wax  spheres  immersed 
in  fluid  (Bagnold,  1954).  However,  because  of  the  geometry  of  Bagnold's 
experiments,  and  because  the  particles  were  neutrally  buoyant,  there  are 
still  many  unanswered  questions  relating  to  granular  shear -flows  of 
geophysical  interest.  We  think  some  of  these  questions  can  be  answered 
by  conducting  a  series  of  basic  laboratory  experiments  designed  to  study 
horizontal  shear  flows  of  negatively  buoyant  particles. 

Accordingly,  an  annular  shear  cell  has  been  constructed  in 
order  to  study  the  relationship  between  the  depth  of  motion  in  a  granula 
bed  and  the  traction  applied  at  the  surface  of  the  bed.  Data  have  been 
collected  for  dry  granular  shear  flows  of  spherical  glass  beads  under 
a  wide  range  of  conditions. 

If  sufficiently  intense  traction  is  applied  at  the  surface  of 
granular  bed,  the  uppermost  grains  dilate  and  shear  relative  to  each 
other.  The  deDth  to  which  motion  occurs  in  the  bed  is  believed  to  be 
governed  by  a  dynamic  yield  criterion  similiar  to  the  well  established 
Coulomb  criterion  for  static  yield.  For  steady  shear  flows  of  cohesion¬ 
less  gra  ties  the  dynamic  yield  criterion  relates  T,  the  applied  shear 


stress,  to  P,  the  normal  stress  at  the  yield  plane,  by  the  relation 
r=P  tan^,  where  tan$r  is  the  dynamic  coefficient  of  internal  friction, 
and  <(>r  is  the  dynamic  angle  of  internal  friction.  When  applied  to  the 
case  of  steady  sediment  transport  the  dynamic  yield  criterion  becomes  a 
relationship  between  the  bedload  and  the  applied  shear. 

Since  commercially  available  viscometric  apparatus  is  totally 
inadequate  for  studying  granular  materials,  a  new  shear  cell  had  to  be 
designed  and  built  (See  Figure  Bl-2).  Shearing  of  the  material  takes 
place  between  two  parallel,  circular,  horizontal  plates.  The  upper  plate 
remains  stationary  while  the  lower  plate  rotates.  The  annular  gap  within 
which  shearing  takes  place  is  4.4  cm  wide  with  a  mean  radius  of  12.4  cm. 

The  shear  cell  was  designed  to  permit  independent  control  of  the  shear 
rate  and  normal  stress.  Shear  rates  and  normal  stress  levels  were  chosen 
to  keep  centrifugal  effects  small.  The  upper  plate  is  free  to  move 
vertically,  thus  allowing  for  d’litation  in  the  granular  material.  The 
outer  wall  is  plexiglass  lined  with  a  removable  sheet  of  clear  polycarbonate. 
This  allows  for  visual  observations  of  the  flow  as  well  as  measurement  of 
the  depth  of  motion.  Thus  the  general  shear  response  for  granular  media, 
as  well  as  the  dynamic  yield  criterion,  can  be  studied  over  a  range  of 
conditions . 

A  series  of  laboratory  experiments  has  been  proposed  to  study 
the  behavior  of  tan<t>r,the  depth  of  motion,  and  the  general  shear  response 
of  granular-fluid  media.  To  date,  data  has  been  collected  for  dry  flows 
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of  spherical  glass  beads  of  densities  p$  »  2.47,2.78,  and  2.97  grams/cm 
and  diameters  1.1  mm  and  1.85  mm.  The  data  verifies  the  general  Bagnoldian 
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concept  that  for  fully  mobilized  shear  flows  at  high  shear  rates,  the 
shear  stress  is  quatratic  in  shear  rate,  as  seen  in  Figure  Bl-3.  This 
behavior  is  quite  different  that  that  of  a  Newtonian  fluid,  for  which 
shear  stress  is  linear  in  shear  rate.  Cursory  examination  of  the  data 
indicates  that  under  many  conditions  there  can  exist  a  shearing  layer  of 
grains  on  top  of  a  nonshearing  layer,  as  expected.  For  these  considtions 
the  parameter  tan<t>r  varies  within  a  narrow  range,  showing  some  dependence 
upon  shear  rate  and  granular  concentration. 

Following  further  analysis  of  this  data  a  similiar  series  of 
experiments  will  be  run  using  water  as  an  interstitial  fluid.  Both  wet 
and  dry  experiments  will  also  be  performed  on  quartz  sand  from  a  natural 
beach. 

We  are  hopeful  that  these  experiments  will  result  in  the 
ability  to  accurately  predict  the  magnitude  of  the  bedload.  Furthermore, 
we  think  the  data  will  prove  quite  valuable  in  the  verification  of  the 
equations  of  motion  for  granular-fluid  media.  Once  these  equations  and 
the  appropriate  boundary  conditions  are  verified,  a  momentum  based  model 
for  bedload  sediment  transport  may  be  formulated. 
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B2.  On-Offshore  Transport  of  Sand  (T.  E.  White  and  D.  L.  Inman) 

The  goal  of  this  study  is  the  prediction  of  on-offshore  transport 
and  the  slopes  of  beach  profiles  from  forcing  functions,  particularly  waves 
and  currents.  This  is  a  field  study  employing  tracer  techniques  and  bed¬ 
load  sensors  (such  as  that  being  developed  in  Task  B1.3)  to  measure  sand 
movement,  and  wave  arrays  and  current  meters  to  measure  the  forcing 
functions. 

Previous  ONR-sponsored  work  (e.g.  Winant  et  al  1975;  Aubrey 
et  al  1980)  focused  on  the  statistical  description  of  beach  profiles,  and 
in  a  qualitative  sense,  related  profile  changes  to  pnysical  forcing 
functions  via  empirical  eigenfunction  analysis.  The  present  study  is 
a  more  direct,  mechanical  approach  to  measuring  the  short-term  (minutes 
to  hours)  bed  response.  Testing  of  state-of-the-art  transport  models  such 
as  that  of  Bailard  and  Inman  (1981)  then  becomes  possible  since  variables 
such  as  tracer  layer  thicknesses  and  currents  that  appear  in  the  models 
can  be  measured  directly. 

An  energetics  model  seeks  an  integral  property,  the  total  trans¬ 
port,  by  considering  the  energy  flux  and  the  rate  of  working  for  the  entire 
sample  grid.  These  energy  fluxes  and  rates  of  working  are,  in  turn,  inte¬ 
gral  properties  that  reflect  the  fluid-granular  particle  dynamics  and 
kinematics  summed  over  the  sampling  grid.  Since  the  property  we  seek, 
i.e.  sediment  transport  rate,  is  an  integrated  property,  conceptually 
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and  methodologically  it  makes  good  sense  to  use  energetic  models. 

Energetics  models  relate  local  sediment  transport  rates,  i  , 
to  local  fluid  velocities,  u  ,  through  coefficients  of  efficiency,  e  , 
and  friction,  c^.  An  energetics  bedload  model  was  developed  by  Bailard 
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and  Inman  (1981)  based  on  a  modification  of  Bagnold's  (1963)  steady 
stream-transport  model,  that  allows  for  varying  strengths  of  longshore 
currents  and  varying  incident  wave  angles.  The  Bailard  and  Inman  model 
predicts  local  bedload  transport  as: 


where  <  >  means  time-averaged,  p  is  the  fluid  density,  c^  is  the 

fluid  drag  coefficient,  is  the  efficiency  of  bedload  transport,  tan$ 
is  the  coefficient  of  granular  friction,  B  is  the  local  bed  slope  in  the 
negative  x-direction,  n^  is  unit  vector  in  the  x-direction  (positive 
onshore)  and  u^  is  the  time-vary*'ng  velocity  vector  composed  of  the 
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steady  current  u  and  the  oscillatory  current  u- 

In  addition  to  bedload,  a  total  load  transport  model  requires  a 
suspended  load  component.  An  energetics  model  for  suspended  load  is  being 
developed  along  the  concepts  of  Bagnold  (1963).  However,  differences 
between  Bagnold's  unidirectional  stream  transport  and  that  under  oscil¬ 
latory  waves  introduces  important  modifications,  making  such  a  model  more 
tentative  than  the  bedload  model  for  which  there  is  now  considerable  body 
of  support  (e.g,  Inman  and  Bowen,  1963;  Inman,  Komar  and  Bowen,  1968; 

Komar  and  Inman,  1970;  Bailard  and  Inman,  1979;  etc.). 

Advecti on-diffusion  models  consider  the  kinematics  of  individual 
particles,  integrated  over  all  space  to  obtain  total  transport.  This 
integration  is  extremely  sensitive  to  local  near-bed  concentrations, 
velocities,  and  the  position  of  the  bed  itself,  both  in  terms  of  moving 
and  static  layers.  Thus  there  are  discrete  properties  which  have  large 
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gradients  near  one  of  the  limits  of  integration,  namely  the  bed  whose 
position  is  not  fully  known.  Since  we  seek  an  integrated  property,  it 
is  a  more  logical  procedure  to  calculate  it  from  other  integrated 
properties,  as  is  done  in  energetics  models.  Further,  the  near  bed 
variabilities  introduce  large  errors  into  the  integration  for  the  total 
sediment  transport.  The  value  of  the  advection-diffusion  model,  beyond 
the  obvious  advantage  of  an  independent  estimate  of  transport,  is  that 
it  provides  much  greater  detail  about  depositional  geometry,  and  hence 
gives  more  insight  into  the  mechanics  of  transport  under  different  surf 
and  near-surf  conditions,  and  into  the  effects  of  bedforms  such  as 
ripples  and  cusps  on  the  transport  processes.  A  common  diffusion- 
advection  model  is  obtained  by  modifying  the  diffusion  equation  to 
include  advective  terms  (e.g.,  Ranee,  1963): 


3M 

at 


where  M  is  the  mass  concentration  of  tracer,  k  ,  k  ,  k  are  the 

x  y  2 

diffusion  coefficients  in  the  two  horizontal  directions  and  vertical 
direction  respectively,  and  U(z)  and  V(z)  are  the  horizontal  trans¬ 
port  velocities  which  are  a  function  of  z. 

It  must  be  recognized  that  the  advection  rates  in  the  above 
relation  vary  with  distance  above  the  bed,  and  cannot  be  integrated 
without  this  knowledge.  Once  the  diffusion  coefficients  and  the  advec- 
tions  are  determined  the  equation  can  be  tested  against  the  distribution 
of  tracer  determined  experimentally. 


In  the  fall  of  1980,  four  experiments  were  performed  outside 
the  breakers  at  Torrey  Fines  Beach  to  obtain  data  with  which  to  test 
these  models.  The  principle  elements  of  this  experiment  were  current 
meters  to  measure  the  fluid  velocity,  core  samples  to  measure  sand 
tracer  distribution,  and  reference  rods  to  measure  net  sand  level  changes 
The  core  samples  were  located  in  a  spatial  grid  of  8-12  samples,  and 
7-10  of  these  sample  sets  were  taken  during  each  experiment.  Each 
spatial  sample  set  can  be  used  to  determine  the  transport  rate  during 
the  time  period  between  sampling  of  that  set  and  the  previous  set.  Fur¬ 
thermore,  each  series  of  samples  taken  at  one  location  but  at  successive 
times  yields  a  single  temporal  estimate  of  transport  (Inman  et  al.,  1980) 
thus  5  sets  of  10  samples  each  yields  15  transport  estimates.  These 
estimates  are  obtained  from  the  relation 

T  =  (ps  -  0)  9  v0! 

where  N„  is  the  "at  rest"  solids  concentration,  is  the  thickness 
o  o 

of  the  layer  of  motion  (obtained  from  the  vertical  tracer  profile  in 
cores),  and  U  is  the  transport  velocity,  determined  from  the  relation 

* 

u-— *r*r~ 

where  M  is  the  mass  concentration  of  tracer,  and  x  is  the  vector 
distance. 

Preliminary  analysis  of  the  data  shows  interesting  differences 
between  transport  in  surfzone-wide  tracer  experiments  (e.g.,  Komar  and 
Inman,  1970;  Inman  et  al.,  1980)  and  those  estimates  obtained  from  local 


experiments  several  meters  outside  the  surfzone.  Outside  the  breakpoint, 
transport  velocity,  U  ,  is  consistently  of  the  order  of  cm  per  minute, 
averaging  about  5  cm  min1  for  onshore  transport  and  2  cm  min1  for  long¬ 
shore  transport.  By  comparison,  surfzone  experiments  show  longshore 
transport  velocities  30  to  60  times  greater. 

Thus,  the  assumption  that  longshore  transport  rates  are  much 
higher  inside  the  surfzone  than  outside  appears  to  be  valid.  But  the 
offshore  transports  show  that  surfzone  transport  is  not  total  transport 
and  the  prediction  of  beach  profile  changes  requires  knowledge  of  the 
offshore  transport  relations.  Estimates  of  transport  rates  also  require 
that  ZQ  ,  the  thickness  of  the  layer  of  motion  be  known.  Various  sta¬ 
tistical  estimators  of  ZQ  ,  which  use  vertical  tracer  concentration 
profiles  as  the  input,  have  been  proposed  (Crickmore,  1967;  Inman  et 
al.,  1980;  Kraus  et  al . ,  1981).  For  complicated  concentration  profiles, 
namely  the  ones  obtained  in  surfzone  experiments,  the  more  accepted 
estimators  vary  among  themselves  by  about  10-20%.  However,  in  the  local 
offshore  experiments  at  Torrey  Pines,  concentration  profiles  generally 
exhibited  abrupt  vertical  cutoff  of  tracer,  thereby  yielding  unambiguous 
estimates  of  the  thickness  of  the  layer  of  motion.  These  thicknesses 
could  often  be  measured  in  the  core  samples  with  a  ruler,  because  of 
high  tracer  concentrations. 

Another  difference  between  the  two  types  of  experiments  is 
revealed  by  the  actual  magnitude  of  the  estimates  of  ZQ  .  Estimates 
from  offshore  experiments  usually  give  thicknesses  of  1-2  cm.  Surfzone 
experiments  sometimes  yield  thicknesses  that  often  range  up  to  7-10  cm. 

We  suggest  that  the  difference  in  magnitude  and  variance  of  ZQ  across 


the  surfzone  is  related  to  the  differences  in  bottom  velocities  under 
plumging  waves,  bores,  and  backwash  vortices  repsectively.  Plunging 
breakers  have  intense  downward  motion  near  the  bottom  that  scours  the 
sand  bed,  as  do  the  intense  vortices  in  the  backwash.  In  contrast  the 
bores  that  traverse  the  surfzone  produce  a  horizontal  velocity  and 
tangential  stress  much  like  a  progressive  wave.  Accordingly  the 'measured 
values  of  ZQ  are  greatest  near  the  breakpoint  and  in  the  swashzone  and 
least  in  the  central  portion  of  the  surfzone. 

The  sand  tracer  data  has  been  analyzed,  yielding  the  transport 
velocities  U  ,  and  the  thickness  of  the  layer  of  motion  ZQ  ,  and  the 
transport  rate  from  their  product.  Once  the  velocity  field  of  the  water 
is  obtained  by  analysis  of  the  current  meter  data  then  all  of  the  neces¬ 
sary  variables  will  be  available  and  transport  model  testing  can  proceed. 

Also,  additional  field  experiments  are  being  conducted  in  con¬ 
junction  with  the  study  of  oscillatory  bursting  described  in  Task  Bl.l. 
Finally,  comparisons  of  the  thickness  of  the  moving  layer  obtained  in 
Task  Bl.l,  in  the  present  Task  B2,  and  with  estimates  based  on  the 
laboratory  study  of  rheology  (Task  B1.2)  will  provide  valuable  control 
limits  as  well  as  considerable  insight  into  the  mechanics  of  fluid- 
granular  interactions. 


B3 .  Coastal  Forms 


Research  has  continued  on  this  task,  which  was  first  outlined  at 
the  ONR  Coastal  Dynamics  Conference,  Panama  City,  Florida.  The  goal  is 
the  development  of  a  classification  criterion  for  coastal  forms  based 
on  the  nature  of  the  fluid  driving  forces  and  the  regime  of  the  relative 
grain  stress  near  the  bed  surface  (Task  3.1).  In  this  model,  the  nature 
of  the  fluid  motion  and  the  resulting  fluid-granular  interactions  deter¬ 
mine  whether  the  bedform  is  of  the  "interactive"  type,  and  free  to  grow 
in  size  as  under  unidirectional  flow  leading  from  ripples  to  dunes;  or 
of  the  "forced"  type  where  the  bedform  geometry  is  constrained  by  length 
scales  in  the  fluid  motion  as  in  the  case  of  oscillatory  waves  where  the 
orbital  displacement  determines  the  wavelength  of  oscillatory  ripples. 

B3.1  Bedforms  Produced  by  Fluid-Granular  Interaction  (D.  L.  Inman  and 

S.  A.  Jenkins) 

Observations  continue  to  verify  the  generality  of  the  genetic  classi¬ 
fication  of  bedforms  into  "interactive"  and  "forced"  forms.  The  classi¬ 
fication  scheme  is  shown  in  Figure  B3-1  and  a  summary  of  the  concepts 
leading  to  the  primary  types  of  bedform  follows. 

Fluid-granular  interactions  generate  two  primary  classes  of  bed- 
forms.  The  first  is  a  class  whose  scales  are  determined  by  strong  feedback 
mechanisms  between  bed  perturbations  and  fluid  eddies  The  scales  of  these 
forms  cascades  upward  from  small  to  larger  features  as  grains  to  ripples 
to  dunes,  somewhat  analogous  to  the  vortex  pairing  in  a  turbulent  shear 
layer.  Since  the  fluid  forming  the  features  is  essentially  unbounded  by 
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Figure  B3-1. Classification  of  bed  forms  in  terms  of  fluid  motion  and  the 

relative  grain  stress.  Numbers  1  through  4  indicate  the  rough¬ 
ness  scales  of  "interactive"  bed  forms  whose  relative  size 
depends  primarily  upon  the  upward  cascade  of  eddy  shedding  per¬ 
turbations  of  unconstrained  fluid  flow  over  an  interactive  bed. 
The  asterisks  *  indicate  "forced"  forms  whose  scale  and  geometry 
are  primarily  controlled  by  the  constraints  imposed  by  the 
wave-like  nature  of  the  fluid  motion;  (2)  and  (3)  following  the 
form  indicate  two  or  three  dimensionality  of  the  form. 
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fluid  constraints  except  at  the  bed  where  multiple  scale  interactions  be¬ 
tween  the  bed  and  the  fluid  occur,  the  bed  features  are  referred  to  as 
"interactive  forms" .  Typical  interactive  forms  include  ripples  and  dunes 
generated  by  the  unidirectional  flow  of  wind  and  water. 

The  scale  and  shape  of  the  second  class  of  bedforms  are  forced  by 
exterior  bounding  constraints  imposed  upon  the  fluid,  such  as  the  wave¬ 
length  of  surface  waves  and  the  orbital  diameters  of  fluid  motions,  and 
are  referred  to  as  "forced  forms".  Feedback  in  forced  forms  is  initially 
absent  and  only  occurs  after  the  bed  feature  becomes  well  formed.  Typical 
forced  forms  include  wave-induced  ripples,  sand  waves  (anti -dunes ) ,  and( 
beach  features  such  as  cusps,  crescentic  bars  and  backwash  ripples. 

Although  forced  forms  have  no  initial  requirement  of  feedback  from  the  bed 
to  the  fluid,  feedback  becomes  progressively  more  intense  as  the  bedform 
grows  in  amplitude.  Eventually  the  bedform  feedback  becomes  the  limiting 
factor  in  determining  the  equal ibri urn  amplitude  of  the  bedform.  The 
uniformity  in  heights  of  oscillatory  vortex  ripples  and  of  beach  cusps 
are  examples  of  such  controlling  feedback. 

The  inertia  accumulated  in  large  scale  bedforms  sometimes  results 
in  their  persisti ng  long  after  the  driving  force  that  generated  them  has 
ceased.  Obvious  examples  are  longshore  bars,  rip  channels,  and  beach  cusps, 
which  because  of  their  great  inertia  may  control  the  fluid  motion  for  some 
time  after  the  wave  regime  that  formed  them  has  changed.  In  these  cases 
the  features  may  be  referred  to  as  "remnant  forced  forms". 

Beach  cusps,  a  case  where  the  amplitude  of  a  forced  form  is  limited 
by  feedback  is  described  in  the  following  section. 
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B3.2  Cusps  (R.  T.  Guza) 

Beach  cusps  are  crescentic  shoreline  features  having  pseudo-regular 
longshore  wavelength  ranging  from  a  few  centimeters  to  as  much  as  several 
hundred  meters.  Cusp  heights  (horn  to  trough  vertical  distance)  vary 
from  a  few  centimeters  to  more  than  2  meters.  In  the  context  of  the  bed- 
form  classification  scheme  discussed  above,  beach  cusps  are  forced  forms. 
The  cusp  wavelength  is  necessarily  equal  to  half  that  of  the  edge  wave 
responsible  for  cusp  growth.  Initial  cusp  growth  is  simply  the  sediment 
response  to  stresses  arising  from  the  longshore  periodic  flow  field  of 
the  edge  wave. 

Progress  has  been  made  in  investigating  the  feedback  from  cusps  to 
edge  waves  which  occurs  when  the  cusp  has  grown  to  large  amplitude.  An 
equilibrium  cusp  amplitude  is  found  such  that  the  tendency  of  the  incident 
waves  to  erase  the  beach  cusps  is  balanced  by  the  edge  waves  trying  to 
increase  the  cusp  amplitude.  The  theory,  which  uses  the  formalism  of  non¬ 
linear  triads,  predicts  that  cusp  heights  depend  on  the  mean  beach  slope 
and  incident  wave  periods  and  amplitudes.  The  model  prediction  can  also 
be  cast  into  a  form  where  the  cusp  amplitude  depends  only  on  the  cusp 
wavelength,  mean  beach  slope,  and  incident  wave  amplitude.  In  the  most 
idealized  case  where  the  presence  of  cusps  does  not  alter  the  mean  beach 
slope  from  its  value  without  cusps,  and  with  strongly  reflected .incident 
waves,  theory  predicts 


<  0.1 


with  6  the  mean  beach  slope  in  radians,  A  the  cusp  wavelength  and 
a^  equal  to  half  the  cusp  heights.  Mo  defines  a  minensionless  number 
which  the  theory  suggests  is  the  most  important  parameter  in  the  problem. 

It  turns  out  to  be  the  ratio  between  the  cusp  steepness  (a^/A)  and  the 
slope  of  the  beach. 

Note  that  the  above  equation  predicts  a  maximum  possible  value  for 
Mo,  so  that  given  the  mean  beach  slope  6  and  cusp  wavelength  A  ,  all 
cusp  amplitudes  for  which  Mq<0.1  are  possible.  The  theory  is  actually 
capable  of  predicting  precise  values  for  Mq,  but  this  requires  detailed 
information  about  the  incident  wave  field  which  is  not  generally  available. 

-To  test  the  theoretical  predictions  observations  of  cusp  amplitude 
ay  were  plotted  versus  as  (Figure  B3-2).  The  data,  which  cover  cusp 
amplitudes  from  a  few  centimeters  to  one  meter,  fall  surprisingly  close  to 
the  maximum  vain  3  predicted  by  theory,  Mq  =  0.1.  Although  the  theory 
makes  many  approximations  and  simplifications,  the  field  data  bears  out 
that  Mq  is  indeed  an  important  nondimensional  number.  Details  of  the 
theory  and  comparisons  to  field  data  are  given  in  Guza  and  Bowen  (1981). 
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B4.  Shore  Processes  of  Northern  Adriatic  Beaches  (m.  Weydert  and  D.  L.  Inman) 

The  long-term  scientific  objective  of  this  cooperative  project 
is  to  perform  a  detailed  study  of  coastal  bathymetry  and  shoreline 
changes  in  a  shallow  water  region  where  wind  waves,  tidal  Kelvin  waves 
and  local  wind-generated  circulation  all  are  important  driving  forces. 

The  study  is  in  a  setting  where  the  historic  shoreline  changes  are  gen¬ 
erally  known  since  800  B.C.  and  well  known  since  the  first  century  A.D. 

Much  of  our  knowledge  of  historical  changes  results  from  the  fact  that 
the  coastal  city  of  Ravenna  was  the  Adriatic  Port  for  the  Roman  fleet 
since  the  first  century  A.D.,  so  that  the  Roman  Legion's  Corps  of 
Engineers  made  a  number  of  harbor  plans  and  conducted  surveys  of  the 
coastline.  Analysis  of  modern  shoreline  changes,  with  special  emphasis 
on  the  sediment  supply  by  rivers  and  the  distribution  of  sediment  along 
the  coast  by  waves  and  currents  will  eventually  yield  a  budget  of  sedi¬ 
ment  for  the  Emilia  Romagna  littoral  cell. 

As  a  first  step  toward  defining  the  offshore  wave  climate  in 
the  northern  Adriatic  region,  a  wave  and  current  array  was  installed  on 
an  AGIP  gas  production  platform  off  the  coast  of  Ravenna,  Italy  in 
July  1980.  Fouling  and  damage  (Figure  B4-1)  of  the  electromagnetic 
current  meter  was  repaired  in  July  1981.  The  available  data  will  be  ana¬ 
lyzed  for  gross  frequency-directional  spectra  and  spectral  moments  using 
methods  discussed  under  Task  Al.  The  wave  data  will  be  refracted  into 
shore  and  used  to  predict  the  strength  and  direction  of  longshore  currents. 

Another  method  for  predicting  the  wave  climate  and  coastal  cur¬ 
rents  is  being  implemented  by  our  Italian  counterparts.  Cavalieri  and 
Rizzoli  (1977)  have  constructed  a  numerical  wind-wave  generation  model 


Figure  B4-1.  Underwater  view  of  fouled  electromagnetic  current  meter  probe 
mounted  to  AGIP  tower  off  Ravenna,  Italy.  Probe  support  is  in 
fact  broken,  but  being  held  in  place  by  large  mussel  growth. 
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used  to  compute  wave  statistics,  longshore  currents  (Figure  B4-Z)  and 
other  variables  in  the  northern  Adriatic  from  meteorological  observa¬ 
tions  on  the  Italian  and  Yugoslavian  coast.  The  model  has  been  used 
to  compute  wave  statistics  furing  the  periods  11  to  18  July  1978  and 
15  September  to  16  October  1978.  Work  is  underway  to  generate  the 
model  predictions  during  the  time  of  actual  wave  measurements  in  1980 
and  1981  so  that  the  two  methods  may  be  compared. 

The  wave  and  current  measurements  and  predictions  will  also 
be  used  to  estimate  the  dissipation  and  scattering  of  wave  energy  in 
this  shallow  sea.  Definitive  results  in  this  area  will  depend  upon 
the  analysis  of  the  data,  currently  under  study. 


Figure  B4-2.  Longshore  current  prediction  from  wind-wave  generation  model.  The  spacing  between  each 
station  number  on  vertical  axis  Is  current  speed  of  1  meter  per  sec.  measured  from  a 
line  drawn  through  the  station  number.  Up  Is  northerly  current.  Note  the  storm  event 
1-4  October  1978. 


C.  COASTAL  ZONE  REMOTE  SENSING 


The  objective  of  the  remote  sensing  task  is  to  use  and  refine  existing 
techniques  to  study  beaches  and  shelves,  rather  than  to  develop  new  remote 
sensing  instrumentation.  The  opportunities  for  wide  spatial  coverage, 
without  great  time  lags,  of  course  makes  remote  sensing  an  attractive  way 
of  gathering  data  which  would  be  prohibitively  expensive  to  obtain  using 
traditional  "ground  truth"  sensors. 

Task  Cl  addresses  the  advances  made  in  SAR  technology  and  suggests 
future  development  work  necessary  for  oceanographically  useful  measurements. 
Tasiv  C2  summarizes  the  data  manipulations  performed  on  satellite  infrared 
images  of  the  California  Current  System  and  the  processing  being  done  to 
compare  these  images  to  current  meter  and  drifter  records. 


Cl.  Remote  Sensing  of  Waves  (SAR)  (S.  S.  Pawka) 

In  the  past  year,  we  have  investigated  the  need  for  future  devel¬ 
opmental  work  on  the  measurement  of  wave  frequency-directional  spectra  with 
synthetic  aperture  radar  (SAR)  systems.  Remote  sensing  techniques  for  wave 
measurement  were  discussed  extensively  at  the  Workshop  on  Wave  Measurement 
Technology,  recently  conducted  by  the  National  Academy  of  Sciences.  The 
Shore  Processes  Laboratory  was  represented  by  Steve  Pawka,  an  invited 
participant.  Our  input  into  the  SAR  discussions  at  the  workshop  drew  heavily 
on  our  past  experience  with  the  comparison  of  L-band  SAR  and  linear  array 
data  sampled  at  Torrey  Pines  Beach,  described  in  detail  in  last  year's 
progress  report. 

Several  conclusions  were  drawn  from  the  discussions  on  SAR  sys¬ 
tems:  1)  Many  advancements  (e.g.  image  enhancement  and  processing  techniques) 
have  been  made  over  the  L-band  technology  used  in  the  West  Coast  Experiment, 

2)  A  majority  of  the  ground  truth  experiments  for  SAR  and  other  radar  systems 
have  been  done  with  low  directional  resolution  surface  sensors  (e.g.  pitch/ 
roll  buoy)  which  do  not  adequately  test  the  high  resolution  of  the  remote 
sensors,  3)  SAR  images  show  promise  of  accurate  wave  directional  measurements 
but  have  an  unknown  transfer  function  which  would  enable  the  estimation 
of  frequency  spectra,  4)  more  ground  truth  verification  of  the  SAR  system 
is  required. 

The  ideal  location  to  stage  a  ground  truth  experiment  is  in  the 
deep  ocean,  where  the  SAR  measurements  are  made.  This  would  eliminate  the 
necessity  to  account  for  local  shelf  refraction  effects  in  relating  SAR 
measurements  to  shallow  water  ground  truth.  However,  no  practical  high 
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resolution  directional  system  has  been  developed  for  use  on  the  surface  in 
deep  water.  Ultimately,  the  most  important  use  of  the  SAR  system  will  be 
in  the  satellite  sampling  of  intense  wave  activity  in  the  deep  ocean. 

Torrey  Pines  Beach  is  on  a  partially  sheltered  coast  (see  Progress  A. 2) 
and  has  relatively  mild  wind  and  wave  conditions.  It  will  not  be  certain 
whether  testing  in  this  location  will  be  applicable  to  deep  ocean  storm 
conditions  until  the  relationship  between  wave  height  and  SAR  image  intensity 
(the  transfer  function)  is  better  understood.  Future  ground  truth  exper¬ 
imentation  at  this  site  should  be  considered  only  after  more  developmental 
work  is  completed  on  th  estimation  of  wave  heights  (frequency  spectra) 
from  SAR  images.  Therefore,  no  SAR/Torrey  Pines  ground  truth  experiments 
are  planned  for  the  immediate  future. 

Further  developmental  work  on  SAR  systems  may  indicate  the  need 
for  ground  truth  experiments  in  mild  environmental  conditions.  Torrey  Pines 
Beach  would  be  a  relatively  desirable  site  for  this  experimentation.  The 
local  shelf  is  narrow  and  relatively  plane.  The  directional  accuracy  of 
array  measurements  and  refraction  calculations  at  this  site  have  been  estab¬ 
lished  at  1-2°  (deep  water  values).  Additionally,  a  variety  of  complicated 
directional  spectrum  forms  are  persistent  in  the  local  wave  climate. 

An  intense  SAR  sampling  pattern  in  the  area  of  the  ground  truth 
sensors  is  recommended.  At  least  10  independent  realizations  of  the  wave 
field  on  each  sample  day  are  required  for  reasonable  comparison  statistics. 
With  a  time  lag  of  roughly  17  minutes  between  independent  samples,  this  would 
require  a  minimum  of  3  hours  of  sampling  at  a  specific  site.  Also,  var¬ 
iable  heading  flight  paths  should  be  employed  to  test  the  directional 
sensitivity  of  the  SAR  images.  The  sampling  should  span  a  period  of  about 
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20-30  days  in  February  or  March  (season  of  high  waves)  to  insure  a  variety 
of  wave  conditions  at  all  wave  frequencies.  It  ’-s  stressed  that  the  utility 
of  wave  measurement  systems  for  coastal  science  or  enginee  ing  requires 
usable  directional  information  in  the  frequency  range  0.05-0.30  Hz. 


Collection  and  preliminary  processing  for  a  data  set  of  infrared 


images  in  the  California  Current  System  from  the  NOAA-6  satellite  is 
approximately  two-thirds  complete.  Sorting  through  data  for  cloud- 
free  areas  has  been  unexpectedly  time-consuming;  however,  recent 
software  improvements  at  the  SRSF  have  decreased  the  time  required 
substantially.  More  than  120  images  have  been  corrected  for  distortion, 
temperature-calibrated  and  mapped  to  a  common  projection.  Sorting  and 
mapping  for  the  entire  data  set,  a  one-year  time  period,  will  be 
completed  in  the  fall  of  1981. 

An  empirical  correction  for  the  effects  of  moisture  in  a  clear 
atmosphere,  developed  by  Dr.  Robert  Bernstein,  will  be  applied  to  the 
images  after  two  filtering  processes  have  been  completed.  The  atmo¬ 
spheric  correction  algorithm  requires  the  difference  of  apparent  tempera¬ 
tures  in  bands  3  and  4  of  the  same  image.  Band  3  (13.7um)  contains 
sensor  noise  equivalent  to  several  tenths  of  a  degree,  while  band  4 
(11  urn)  is  virtually  noise-free.  (See  Figure  C2-1A).  The  difference 
image  contains  information  about  the  atmosphere  and  the  sensor  noise 
from  band  3.  (See  Figure  C2-1B).  Data  which  are  contaminated  by 
clouds  are  apparent  in  this  difference  image  as  light  areas.  Note  that 
the  land,  the  Baja  California  coastline  (right  side  of  image),  and 
Guadalupe  Island  (in  the  center  of  the  image)  are  also  visible  in  the 
difference. 


Two  filters  have  been  developed  which  use  this  difference  image. 
The  first  filter  isolates  the  noise  from  the  AVHRR  (See  Figure  C2.1C). 

The  noise  is  subtracted  from  the  original  band  3  image  to  produce  the 
smoothed  version,  shown  in  Figure  C2-1D.  Removing  the  noise  in  this 
way  does  not  smooth  ocean  features,  although  it  may  blur  somewhat  the 
cloud  edges.  The  second  filter  flags  the  cloudy  areas.  This  filter 
is  an  adaptation  of  the  method  used  by  Coakley  and  Bretherton  (in  press) 
to  determine  cloud-cover  statistics  from  the  same  data.  Because  the 
atmospheric  correction  is  valid  only  for  a  clear  atmosphere,  data 
which  are  contaminated  by  clouds  must  be  excluded  from  the  statistical 
analysis.  All  of  the  processing  procedures  are  nearly  automatic  because 
of  the  large  amount  of  data.  These  filtering  processes  will  be  described 
in  detail  in  a  forthcoming  paper. 

Satellite  data  for  the  Coastal  Ocean  Dynamics  Experiment  (spon¬ 
sored  by  NSF)  area  have  been  collected  since  April  1981.  Preliminary 
analyses  of  doppler  velocities  and  drifter  tracks  indicate  unexpectedly 
energetic  small  scale  motions.  Satellite  data  will  be  useful  in  inter¬ 
preting  these  surface  velocity  measurements .  Images  for  the  CODE  area 
from  late  March  through  the  end  of  June  1981  have  been  processed  through 


the  mapping  step. 
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